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ABSTRACT 

Context. Transition disks are believed to be the final stages of protoplanetary disks, during which a forming planetary system or pho- 
toevaporation processes open a gap in the inner disk, drastically changing the disk structure. From theoretical arguments it is expected 
that dust growth, fragmentation and radial drift are strongly influenced by gas disk structure, and pressure bumps in disks have been 
suggested as key features that may allow grains to converge and grow efficiently. 

Aims. We want to study how the presence of a large planet in a disk influences the growth and radial distribution of dust grains, and 
how observable properties are linked to the mass of the planet. 

Methods. We combine two-dimensional hydrodynamical disk simulations of disk-planet interactions with state-of-the-art coagula- 
tion/fragmentation models to simulate the evolution of dust in a disk which has a gap created by a massive planet. We compute images 
at different wavelengths and illustrate our results using the example of the transition disk LkCal5. 

Results. The gap opened by a planet and the long-range interaction between the planet and the outer disk create a single large pres- 
sure bump outside the planetary orbit. Millimeter-sized particles form and accumulate at the pressure maximum and naturally produce 
ring-shaped sub-millimeter emission that is long-lived because radial drift no longer depletes the large grain population of the disk. 
For large planet masses around 9 Mj„p, the pressure maximum and, therefore, the ring of millimeter particles is located at distances 
that can be more than twice the star-planet separation, creating a large spatial separation between the gas inner edge of the outer disk 
and the peak millimeter emission. Smaller grains do get closer to the gap and we predict how the surface brightness varies at different 
wavelengths. 

Key words, accretion, accretion disk - circumstellar matter -stars: premain-sequence-protoplanetary disk-planet formation - stars: 
individual (LkCal5) 



1. Introduction 

Circumstellar disks are the birthsites of planets. The physical 
conditions and the evolution of these disks control the planet 
formation mechanisms. An impoitant goal is to provide theoret- 
ical models and observational constraints to understand the var- 
ious stages in the evolution of gas and dust in the disk. Decrease 
of mass accretion rate (Sicilia-Aguilar et al. 2010; Fedele et al. 
2010), and near-infrared excess with time (Hernandez et al. 
2007; Andrews et al. 2011b) indicate that disks have a range 
of lifetimes from 1 to 10 Myr. 

With the advent of powerful infrared space telescopes such 
as Spitzer, a new class of objects has been identified, called the 
transition disks (e.g., Espaillat et al. 2010). Their spectral energy 
distribution (SED) and direct sub-millimeter (mm) imaging sug- 
gest that warm dust in the inner disk is strongly depleted com- 
pared to the outer disk. The small number of transition disks 
(Muzerolle et al. 2010) suggests an inside-out evolution that 
occurs rapidly. Various mechanisms have been proposed so far 
to explain the inner disk clearing: photoevaporation winds (e.g. 
Owen et al. 201 1), grain growth (Klahr & Henning 1997; Tanaka 
et al. 2005) and dynamical interactions with companions (Lin & 
Papaloizou 1979). Transition disks are therefore excellent labo- 
ratories for planet formation models. 



The clearing of a gap by a companion or planet, from a sim- 
plistic point of view, depends on the competition between the 
viscous torque from the disk and gravitational torques from the 
planet. For a laminar disk, a 1 Mjup planet can clear a gap or 
hole. The recent discovery of a companion inside a massive disk 
in T Cha (Huelamo et al. 2011) supports the scenario of a dy- 
namical clearing, at least for this object. However, models show 
that a single planet seems unlikely to be capable of creating the 
observed large holes, which require multiple systems (Zhu et al. 
201 1; Dodson-Robinson & Salyk 201 1). 

Interestingly, for most of these transition disks the inner cav- 
ity is not empty. They still present relatively high accretion rates 
(~ lO'^M^yr^; see e.g.,Calvet et al. (2005); Espaillat et al. 
(2007)) which implies that the inner cavity is not completely 
empty and that some gas flows through the gap. To allow mass 
flowing, a limit for the planet mass can be inferred, depending 
on the disk viscosity (Lubow & D'Angelo 2006). In addition to 
the gas, some transition disks also present a strong near-infrared 
excess, indicating the presence of dust close to the star. Rice 
et al. (2006) studied the filtration of dust in the gap, considering 
a fixed size for the dust particles, and concluded that increasing 
the planet mass from e.g. 0.5 Mjup to 5 Mjup, the maximum par- 
ticle size that sweeps into the gap decreases from ~ 10 pm to a 
few tenths of a micrometer. 
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Transition disks are potentially interesting laboratories to 
study processes related to the impact of planet formation on 
the disk. One of the most stubborn problems in planet forma- 
tion is the so-called "meter-size barrier". A one meter size ob- 
ject at 1 AU drifts towards the central star in timescales shorter 
than the growth timescales, impeding it to grow (see e.g., Brauer 
et al. 2008; Birnstiel et al. 2010a). In addition, high relative ve- 
locities lead to numerous fragmentation collisions converting a 
large object into smaller dust particles. The same physical pro- 
cess happens to the millimeter-size particles that are observed in 
the outer regions of the disk (e.g., Wilner et al. 2000; Ricci et al. 
2010; Guilloteau et al. 2011). One possible solution to prevent 
the rapid inward drift and trap dust particles, is to consider pres- 
sure bumps (Klahr & Henning 1997; Fromang & Nelson 2005; 
Brauer et al. 2008; Johansen et al. 2009; Pinilla et al. 2012). A 
long-lived positive pressure gradient can lead dust particles to 
move outwards, causing an accumulation of dust at the location 
of the pressure maximum. A large pressure bump is expected in 
protoplanetary disks as a consequence of the presence of a mas- 
sive planet in a disk. In fact, when a planet carves a gap in a disk, 
the gas surface density shows a significant depletion, resulting 
in a large pressure bump at the outer edge of the gap. The dust 
material is trapped and piles up in this local pressure maximum 
where the gas motion is exactly Keplerian, and as a result there 
is no frictional drag between the gas and the particles. For that 
reason, not only do the particles not drift anymore, they also do 
not experience the potentially damaging high-velocity collisions 
due to relative radial and azimuthal drift. Under those circum- 
stances, growth to larger-than-usual sizes is expected, possibly 
even a breakthrough that leads to overcome the growth barrier 
However, if turbulence is still strong enough, particles may frag- 
ment due to their relative turbulent velocities. This scenario is 
only possible if we assume the presence of a planet, formed by 
another mechanism than dust agglomeration. 

In this paper, our goal is to test the idea that the outer edge 
of a planetary gap is a particle trap. We consider the dust evo- 
lution in a disk, where the gas density profile is determined by 
its interaction with a massive planet. We then explore the case 
of LkCa 15, a transition disk (Mulders et al. 2010). The disk has 
been intensively observed at millimeter wavelengths (Pietu et al. 
2006; Andrews et al. 201 lb) with a maximum angular resolution 
corresponding to 28 AU (Isella et al. 2012). The dust continuum 
images show a ring-like structure from ~ 42 to ~120 AU, which 
is best fitted by a flat surface density profile. In addition, a ~ 6 - 
15 Mjup planet in circular orbit and located at 15.7 ± 2.1 AU, 
was claimed (Kraus & Ireland 2011), but not confirmed yet. 

We describe the numerical simulations of planet-disk inter- 
actions that are considered for this work, as well as the coag- 
ulation/fragmentation model of the dust evolution in Sect. 2. 
Section 3 describes the results of the numerical simulations. In 
addition, we present observation predictions in Sect. 4. Section 5 
is a discussion of the main results of our work and Sect. 6 is a 
summary. 



2. Disk-planet interactions and dust evolution 

For the disk-planet interaction process, we use the two- 
dimensional hydrodynamical code FARGO (Masset 2000), that 
solves the Navier-Stokes and continuity equations. For the dust 
evolution, we use the code described in Birnstiel et al. (2010a), 
that computes the coagulation and fragmentation of dust grains, 
including turbulence, gas drag and radial drift. 



Table 1. FARGO parameters 



Parameter 



Values 



Eo[M*/r?] (flared disk) 
I^olMJrj] (wedge disk) 

a 



rmlrp] 

'"out['' p] 

rir X 

6[h{r)] 
Boundary 
conditions 



1.0 
1.26x10-3 
7.10x10-'' 
{1.0,9.01 
llO-\ lO-\ IQ--] 
{-1,-0.51 

0.1 

7.0 

[512 X 1024] 
{0.6,0.81 
open 



2.1. Codes 

FARGO uses finite differences and a fixed grid in cylindrical 
coordinates (r, 0).lts advantage lies in the calculation of the az- 
imuthal advection. It allows each annulus of cells to rotate at its 
local Keplerian velocity and at the end of each timescale, the 
results are stitched together This technique speeds up the cor- 
responding calculations. FARGO has been used for numerous 
studies of disk evolution (e.g., Crida et al. 2006; Zhu et al. 201 1; 
Regaly et al. 2011). We selected open inner and outer bound- 
ary conditions such that the material can leave the grid when the 
flow velocity at the inner or outer radial cells is pointing out- 
wards. For this work, we only focus on the gas density profile at 
the outer edge of the open gap. 

For the dust evolution, we consider that grains grow, crater 
and fragment due to radial drift, turbulent mixing and gas drag. 
Brownian motion, turbulence, vertical settling, radial and az- 
imuthal drifts are taken into account to calculate the relative 
velocities of dust particles. For micron-sized particles, the rela- 
tive velocities mainly result from Brownian motion and settling. 
These grains grow by coagulation as a result of van der Waals in- 
teractions. As they grow, they start to decouple from the gas, and 
turbulent motion as well as radial drift become the main sources 
for relative velocities. 

When dust grains reach sizes with high enough velocities, 
they encounter destructive collisions. Particles above the frag- 
mentation velocity always cause fragmentation if the impactor 
mass is within one order of magnitude of the target, otherwise 
they cause erosion. The transition from the 100% sticking prob- 
ability to 100% fragmentation/erosion probability is linearly in- 
creasing between 0.2 x Vf to v/, as in Birnstiel et al. (2010a). 

The fragmentation velocities of the dust grains can be 
estimated based on laboratory experiments and theoretical work 
of collisions for silicates and ices (e .g., Blum & Wurm 2008; 
Schafer et al. 2007; Wada et al. 2009). For silicates, they are of 
the order of few m s"' (e .g., Blum & Wurm 2008), and increase 
with the presence of ices (Schafer et al. 2007; Wada et al. 2009). 

To characterize how well the dust particles couple to the gas, 
we refer to the Stokes number St, which is defined as the ratio 
between the turn-over time of the largest eddy (I/O) and the 
stopping time of the particle due to the friction with the gas. In 
the presence of a drag force, big particles are not affected by the 
drag and they move with their own Keplerian velocity, implying 
that St » 1. When St 1, the particles are small enough to 
be well coupled to the gas and move with the gas. 

In the Epstein regime, where the ratio between the mean free 
path of the gas molecules A^fp and the sizes of the particles a is 
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Table 2. Dust model parameters 



Parameter 


Values 




1.0 


Afd,sk[Mj„p] 


55.0 


a 


{lO-\ IQ-'-] 




2.0 


'■ou.[AU] 


140.0 


R. [Ro] 


1.7 


[K] 


4730 


p, [g/cm^] 


1.2 


Vf [cm/s] 


1000 



/Imfp/a ^ 4/9, the Stokes number is at the midplane given by 
Birnstiel et al. (2010a), 



(1) 



where p, is the intrinsic density of the dust particles, and is 
the gas surface density given by 'Lg - p{r, z)dz, where p(r, z) 
is the total gas density. 

Due to the sub-Keplerian velocity of the gas, particles that 
are large enough and do not move with the gas velocity experi- 
ence a head wind, leading them to lose angular momentum and 
move inwards. The radial component of the dust velocity is the 
sum of a term corresponding to the frictional drag force (which 
opposes the motion of the dust grain) and of a term correspond- 
ing to the drift velocity, Mdrfft, with respect to the gas, which di- 
rectly depends on the radial pressure gradient drP- It is given 
by: 



^r.dust ~ 



•*r,gas 



1 + St^ st-1 + St pi^ ' 



(2) 



where Mi gas is the gas radial velocity. The first term is the domi- 
nant contribution for St <k 1 . Equation 2 assumes low values for 
the dust-to-gas ratio. The drift velocity is usually negative, un- 
less the pressure gradient is positive at a given radius. Particles 
with St = 1 experience the highest possible radial drift. 

When particles reach certain sizes due to coagulation, they 
have relative velocities that may be high enough to cause frag- 
mentation. For particles with St ~ 1, the maximum turbulent 
relative velocity is given by (Ormel & Cuzzi 2007) 



(3) 



where a is the turbulence parameter (Shakura & Sunyaev 1973). 

For dust particles with a/2 < St < 1 relative velocities 
are dominated by turbulent motion and radial drift. However, if 
the radial drift is reduced due to a positive pressure gradient, an 
approximation of the maximum size of the particles is obtained 
when the turbulent relative velocities are as high as the fragmen- 
tation velocity v/, therefore: 



42„ 



7 



(4) 



The dust grain distribution n(r, z, a) depends on coagulation and 
fragmentation collisions and is described by the Smoluchowski 
equation (see Birnstiel et al. 2010a, Eq. 35 and Eq. 36). The 
vertically integrated dust surface density distribution per loga- 
rithmic bin is 



cr(r, a) 



xJ-a. 



n(r, z,a) ■ m - adz 



(5) 



where m is the mass of a single particle of size a, hence the total 
dust surface density is: 



Jo 



cr(r,a)dlna. 



(6) 



For the coagulation process, the midplane is the most impor- 
tant region, therefore the vertically integrated dust surface den- 
sity distribution is a good approximation for describing the dust 
grain distribution. The advection-difFusion differential equation 
that describes the evolution of the dust surface density 2^ for a 
single dust size, can be written in cylindrical coordinates as: 



dt 



\_d_ 
r dr 



ld_ 
r dr 



r^gDd 



hi 



0, 



(7) 



where Dd is the dust diffusivity and is equal to 
Dd - v/(l + St^), if the gas diffusivity is taken to be the 
turbulent gas viscosity v - a h. Equation 7 is solved for each 
size using the flux-conserving donor-cell scheme (see Birnstiel 
et al. 2010a, Appendix A). 

For a massive planet on a fixed orbit, the gap opening pro- 
cess reaches a quasi-steady state in timescales that are shorter 
(< 1000 orbits) than the million-years timescales of the entire 
process of dust evolution. For this reason, we use the gas surface 
density after one thousand orbits (~ 10"^ Myr) and take it as an 
input for the dust evolution models. Therefore, we assume that 
the gas surface density remains stationary from one thousand or- 
bits to several million years. 

We consider that after the giant planet is formed, the disk 
remains dust-rich and we study how the remaining dust evolves 
after the gap is open and has reached a quasi-steady state. We 
consider that the gas does not have any feedback from the dust 
and stays static during the dust evolution process since the dust 
density is always less than the gas density. In addition, since 
the orbital timescales are much shorter than the dust evolution 
timescales and there are no strong asymmetries produced by the 
planet, we azimuthally average the surface density for our mod- 
els. With these assumptions, the dust evolution at the outer edge 
of the gap is computed in an accurate way. 

2.2. Set-up 

For our simulations, we assume a locally isothermal disk and 
two different radial density profiles 1, cc r^. When the disk 
scale increases with radius as h/r cc r'^'* (flared disk) as sug- 
gested by modelling of T Tauri disks (D'Alessio et al. 2001), 
the temperature scales as T cc r^'^^, and /3 - -1. When h/r is 
constant (wedge disk), the temperature scales as T oc r"' and 
/? = -0.5. The a turbulence parameter is considered to have val- 
ues between lO"** to 10"^, in agreement with magnetorotational 
instability (MRI) of active disks (Johansen & Klahr 2005). The 
numerical simulations are performed with - 512 radial log- 
arithmically, and - 1024 azimuthal uniformly distributed 
grid cells. 

All parameters such as the stellar properties, kinematic vis- 
cosity, disk mass, the scale height and the gas density profile of 
the disk are taken to be the same between the hydrodynamical 
gas evolution and the dust coagulation simulations. 

For FARGO simulations, we use normalized units such that 
G - Mi, +Mp - 1 and the location of the planet is at r = 1 . The 
simulations are performed from Tin = 0.1 to rout = 7.0, such 
that the planet is at rp - 20 AU for a grid from 2 AU to 140 AU. 
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20 40 50 80 100 120 140 20 40 60 80 100 120 140 



r (AU) r (AU) 

Fig. 1. Azimuthally averaged surface density after 1000 orbits and three different values of a-turbulence for a gap created by a planet 
located in a fixed orbit at 20 AU of 1 Mjup in a flared disk (top-left), 9 Mjup in a flared disk (top-right), 1 Mjup in a wedge disk 
(bottom-left) and 9 Mjup in a wedge disk (bottom-right). The peak of the gas surface density at the position of the planet is due to 
the strong gravitational effect of the massive planet on the gas of the disk. 



The initial gas surface density Zq at the position of the planet is 
fixed S = 'Zo{r/rp'f, such that it is around ~ 1-1.5 the minimum 
mass solar nebula (Hayashi 1981), and the disk mass is around 
~0.055 M*. We consider a planet at a fixed orbit, with two differ- 
ent masses, 1 Mjup and 9 Mjup. The planet is introduced slowly 
into the smooth disk to avoid numerical issues. Planetary accre- 
tion is not taken into account since this process is still debated. 
We use M* - 1 Mq, and consequently, Mdisk = 55 Mjup. 
Finally, the gravitational effect of the planet is smoothed out, 
such that the gravitational potential (p is softened over distances 
comparable to the disk scale height: 



GMp 



(r2 + e2)j 



(8) 



where e is taken to be 0.6 h. The main FARGO parameters are 
summarized in Table 1 . 

For the dust evolution, we consider that all particles are ini- 
tially 1 fim large, and a dust-to-gas ratio is 0.01, with an intrinsic 
volume density of p,. = 1.2 gem"-' and a fragmentation velocity 
ofvy = 1000cm s"'. Two additional stellar parameters are taken 
into account, the stellar radius which is taken to be 7?* - 1 .7 Rq 
and the star effective temperature T^, = 4730 K, as typical values 
for T Tauri stars. The parameters of the dust evolution simula- 
tions are in Table 2. 



3. Results 

3. 1 . Gas density profile 

Figure 1 shows the azimuthally averaged gas surface density af- 
ter 1000 orbits in four different cases: a 1 Mjup and a 9 Mj^p 
planet on a fixed orbit at 20 AU in a flared and wedge disk re- 
spectively. In each case, three different values for a turbulence 
parameter are considered. 

With a 1 Mjup planet and a high value of a = 10"^, the 
viscous torque on the gas exceeds the gravitational torque by the 
planet, resulting in a much less pronounced gap than in the cases 
of a = 10"^ and a = 10 The gap is replenished easily 
because of the high viscosity of the disk as already predicted 
by different authors (e.g de Val-Borro et al. 2006; Crida et al. 
2006). However, for all a values, the width of the gap is roughly 
the same. It is independent of the viscosity and it depends basi- 
cally on the mass and location of the planet. We notice that with 
low viscosity, the density waves that are produced by the planet- 
disk interaction at the gap edges, become more pronounced but 
they eventually damp with radius. Indeed, with a 1 Mjup planet, 
these waves disappear beyond 40 AU and the disk surface den- 
sity profile follows the initial power law. In this case, the main 
difference between the two disk geometries that we consider, is 
that for each value of a, the gap is slightly deeper for a wedge 
than for a flared disk. 

With a 9 Mjup planet, the gravitational torque dominates the 
viscous torque and the depth of the gap is almost independent 
of the values that we consider for a. After 1000 orbits the sur- 
face density reaches almost the same values near the planet and 
the width of the gap slightly changes for the different values 
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r ( AU ) r ( AU ) 

Fig. 2. Radial pressure gradient after 1000 orbits and three different values of a turbulence for a gap created by a planet located in a 
fixed orbit at 20 AU of 1 Mjup and 9 Mjup in a flared disk. 



o 



Q. 



^ 10 



Flared Disk 
1000 orbits 











•— • 
























25 30 35 40 45 50 



r ( AU ) 

Fig. 3. Effect of the planet mass on the radial pressure gradient 
after ~ 1000 orbits in a flared disk, a = 10"^ and five different 
values of the planet mass. The planet is located at 20 AU. The 
vertical lines indicate the locations of the pressure maximum. 



of turbulence, when the disk is less turbulent the gap is mod- 
erately wider. The density perturbation waves are more evident 
in this situation than for 1 Mjup, however, for a - 10^^ - 10"^, 
they again damp with radius, and at around 60 AU the gas sur- 
face density follows the initial power law. After 1000 orbits and 
a - 10 the density perturbation due to the planet is still 
propagating, for this reason the wiggles in the density profile are 
stronger than in the cases of a - 10"^ - 10"^. These wiggles 
are still present for r > 100 AU and they change slightly with 
the geometry of the disk. 

In all four cases, it is possible to distinguish a peak in the sur- 
face density at the position of the planet, due to circumplanetary 
material. This is a natural consequence of the planet insertion in 
a smooth disk, and because the accretion onto the planet is not 
considered, the peak remains unaltered during the simulations. 




3.2. Pressure gradient 

As explained before, the rapid inward drift of the particles stops 
if the gradient of pressure is positive and the particles pile up 
in the pressure maxima. Depending on the strength of the pres- 
sure gradient, the particles can be trapped or not in the pressure 
bump (Pinilla et al. 2012). Figure 2 shows the radial pressure 
gradient considering the azimuthally averaged gas surface den- 
sity after one thousand orbits, in the case of a flared disk and 
with an isothermal equation of state. The pressure is given by 
P(r) = p(r) c^(r), where the sound speed is defined as 



(9) 



where kg is the Boltzmann constant, nij, is the proton mass 
and ju is the mean molecular mass, which is fi - 2.3 in proton 
mass units. If the temperature is considered as a smooth de- 
creasing function with radius, density variations induce pressure 
inhomogeneities. We can see that for a 1 Mjup planet, there is a 
noticeable diff'erence in the amplitude of the pressure gradients 
for a = 10"^ and a - 10^'^ - 10"^. The pressure gradient 
is positive from ~22 AU to ~30 AU, and since Mdnft d,-P 
(Eq. 2), the dust particles have positive drift velocities and move 
outwards. In all cases, the zero point of the pressure gradient is 
located at ~ 10 AU (with a slight variation with viscosity) from 
the planet orbit and does not change with the geometry of the 
disk, since the density profile behaves almost the same for flared 
and wedge disks. 

For a 9 Mjup planet orbiting at 20 AU, the pressure gradi- 
ent on average behaves similarly for each of the a values. The 
wiggly profile of the gradient is due to the density waves that 
are produced by the interaction between the planet and the disk. 
For this planet mass and all values of a, the pressure gradient is 
positive right after the planet location and remains positive until 
around 48 AU, and then becomes again negative for both cases 
of flared and wedge disks. The influence of the planet mass on 
the range where the pressure gradient is positive is remarkable. 
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Flared Disk, Mp =1.0M,„„ & q-10 - 




r ( AU ) 



r ( AU ) 



r ( AU ) 



Fig. 4. Vertically integrated dust density distribution at different times of evolution. Case of a surface density with a gap created 
by 1 Mjup located in a fixed orbit at 20 AU in a flared disk and a = 10"^ (top panel), a - 10^^ (middle panel). Case with no 
planet and a = 10"^ (bottom panel). The solid white line represents the particle size corresponding to St = 1, while the blue line 
corresponds to the fragmentation barrier according to Eq. 4 (Brauer et al. 2008; Birnstiel et al. 2010a). The dashed line corresponds 
to the size for which particles drift faster than those dragged by the gas. Particles above the dashed line are perfectly trapped. 



In both mass regimes, it is clear that the location of the pres- 
sure maximum differs from the outer edge radius of the gap, 
where the gas density again starts to increase. Therefore, if the 
pressure gradient is high enough to moderate the rapid inward 
drift, dust particles are trapped further away from the gap outer 
edge in the gaseous disk. This implies that the gap width for the 
gas is smaller than the one for the dust. 



Figure 3 summarizes these findings, and shows 
pressure gradient curves for additional planet masses, 
Mp = {1,3,5,9, 15)Mjup with orbital radius rp = 20 AU. 
The pressure gradient has its zero point at around 
{29.5, 33.5, 43.5, 48.5, 54.0} AU respectively. 

As a comparison, we use the Hill radius rfj, defined as 
fH - rp {Mp/3Mi,y^^ . In spite of the fact that the outer gap 
edge in the gas is at most 5 rn (see e.g. Dodson-Robinson & 
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Flared Disk, Mp =9.0M,„„ & a = 10"'' 
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Fig. 5. Vertically integrated dust density distribution at different times of evolution for a surface density with a gap created by 9 Mjup 
located in a fixed orbit at 20 AU in a flared disk and a - 10 Lines are as in Fig. 4. 



t=l Myr 




■ - Initial condition 

— 1M,„, 

— 9M,„, 

— 15Af,„„ 



20 40 60 80 100 120 

r ( AU ) 

Fig. 6. Dust-to-gas ratio for 1, 9 and 15 Mjup planet, after 1 Myr 
of dust evolution 



Salyk 2011), we find that the location where dust accumulates 
(i.e., the location of the maximum of the pressure) is at ~ 7 r// for 
1 and 3 Mjup planets, and ~ 10 rn for planets with Mp > 5 Mjup. 
This result will be discussed in more detail in Sect. 5. 

Note that for the FARGO simulations in the case of a 15 Mjup 
planet, we insert the planet in the disk after more orbits than in 
the other cases studied in Sect. 3. 1 , we use non-reflecting bound- 
ary conditions and a smoothing parameter of e = O.S/i (see Eq. 
8) to smooth out the strong density waves produced in this sce- 
nario. 

3.3. Dust size evolution 

For the study of the density distribution of dust particles, we 
again focus on the flared disks, since this disk geometry is in 
good agreement with observations of T Tauri disks (D'Alessio 
et al. 2001). For this section, the surface density peak at the 
planet location was smoothed out to avoid artificial effects at this 
location. 

Top and middle panels of Fig. 4 show the vertically inte- 
grated dust density distribution at different evolution times for 
a gas surface density with a gap created by a 1 Mjup planet lo- 



cated in a fixed orbit at 20 AU, and two different values of a 
turbulence. As explained before the timescales of the gap open- 
ing process are much shorter than the whole dust evolution pro- 
cess, which allows us to take the gas surface density after ~ 1000 
orbits and consider that it remains quasi-stable during the dust 
evolution. The solid white line has the shape of the gas surface 
density profile and is calculated when the particle size coiTe- 
sponds to a Stokes number of unity (Eq. 1). The blue line is the 
maximum size that particles can reach before they have veloc- 
ities higher than the fragmentation limit. It is important to note 
that the maximum size of the particles is inversely proportional 
to the turbulence parameter (Eq. 4), for this reason the maximum 
possible grain size (blue line) differs by one order of magnitude 
for each a turbulent parameter considered. 

The dashed line provides an approximation of the minimum 
size of the particles that are trapped in the pressure maximum. 
Particles with sizes over the dashed line are trapped and the ones 
below the line are dragged by the gas. This condition is found 
when the radial component of the dust velocity (Eq. 2) is posi- 
tive, therefore: 



St > -^pQ, 

drP 



(10) 



which in terms of the particle size and pressure gradient can 
be written as: 



6aS„ 



^critical — 



P Mid log P/d log r)\ 



3 c/log2( 
2 d log r 



(11) 



3.3.1. Effect of turbulence 

We notice a significant influence of the trapping of dust parti- 
cles on the strength of the pressure gradient. Micron-size par- 
ticles (St <s 1), which are easily mixed due to turbulence, are 
more difficult to trap, they go through the gap and finally drift 
towards the star. Dust grains with larger Stokes number can per- 
fectly be trapped because they feel a tail-wind due to the fact that 
the gas moves with super-Keplerian velocity where the pressure 
gradient is positive, resulting in an outwards movement for those 
particles. 

For a = 10"^, the pressure gradient is less steep, the turbu- 
lent mixing and drag along the accretion flow wins over the trap- 
ping, and the dust gets lost onto the star. Even when the pressure 
gradient is positive in a range (~ 22-30) AU, it is not sufficient to 
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Fig. 7. Left: 1.3 mm map in the 1 Mjup case, convolved with a beam of 0.1" x 0.1". Units are in Jy/beam. Right: same for 9 Mjup. 



trap the particles in the pressure maximum. Indeed, after 3 Myr 
there is almost no dust remaining in the whole disk. 

For a = 10"^ turbulence mixing is less strong and in the 
range of (~ 22-30) AU, the pressure gradient is greater (Fig. 2- 
left panel), and is high enough to stop the rapid inward drift of 
the particles. The particles are quickly trapped at the location of 
the pressure maximum (~ 30 AU) and stay there after several 
Myrs. The gap created in this case filters out all particles larger 
than a few fim and the particles that are inside the gap eventually 
drift towards the star The gap formed by 1 Mjup mass planet with 
a = 10"^ produces a pressure gradient high enough to create a 
concentration of dust particles at 10 AU from the location of the 
planet and the particles reach sizes of a few cm. 

As we note in Fig. 2- left panel, for 1 Mjup and a - 10 the 
pressure gradient is slightly higher than in the case of a = 10"^, 
therefore the dust particles should be easily trapped in this case 
reaching even larger sizes (Eq. 4). For comparison, we notice 
in the bottom panel of Fig. 4, that in the case without planet 
and with the same disk and stellar parameters, the dust instead 
of being trapped, grows, fragments, mixes and drifts to the star 
without forming any mm-size pebbles in the outer regions of the 
diskr > 20 AU. 

The necessary size that particles have to reach to be trapped 
(Eq. 1 1) is higher than the maximum sizes reached due to rela- 
tive turbulent motion (Eq. 2). The influence of turbulence for the 
dust dynamics is remarkable, if turbulence is not considered, par- 
ticles grow enough to be decoupled from the gas and they would 
be trapped because of the positive pressure gradient regardless 
of its strength. However, we conclude that the high turbulence 
causes particles to go through the gap and drift finally to the star, 
due to the fact that they reach high relative turbulent velocities 
which favors fragmentation, so the mm-size particles become 
again micron-size grains which are not possible to trap with the 
strength of the pressure gradient produced with or = 10"^. 

3.3.2. Effect of planet mass 

In the case of a 9 Mjup planet, the pressure gradient is roughly 
similar for each of the values of the turbulence parameter a (see 
Fig. 2-right panel), hence for the dust density evolution we con- 
centrate on the intermediate value of the viscosity, a = 10"^. 
Figure 5 illustrates the vertically integrated dust density distribu- 



tion for the case of 9 Mjup at different times of evolution. Similar 
to the case of 1 Mjup planet and a = 10"^, the dust particles 
grow and due to the positive pressure gradient between ~ 22 and 
~ 48 AU (Fig. 2-right panel), the rapid inward drift is counter- 
acted by the positive local pressure gradient and the particles are 
retained in a broad bump of ~ 20 AU width. The dust particles 
in the outer regions of the disk grow until mm sizes and remain 
there even after a few Myrs. It is important to mention that even 
when the pressure gradient reaches a similar value to the case 
of 1 Mjup and a = 10"^, for which there is not particle trapping 
at all, the turbulent relative velocities are lower and the range of 
positive gradient is wider, enabling dust retention (Pinilla et al. 
2012). 

Comparing the dust evolution with a = lO^^* for a 1 Mjup 
(Fig. 4-bottom panel) and a 9 Mjup (Fig. 5), we can conclude 
that the planet mass directly influences the position where the 
dust can be retained and that the dust bump is wider for the more 
massive planet. In addition, the pressure gradient created by the 
different planet masses also affects the size for which particles 
are trapped. For the same turbulence parameter a = 10"^ and 
1 Mjup, the maximum flciiticai is equal to ~ 0.70 cm, while for 
9 Mjup the maximum Ocriticai is ~ 0.47 cm. 



3.3.3. Dust to gas ratio 

It is important to note that for the dust evolution models, we 
assume low values for the dust-to-gas ratio and we do not con- 
sider to have any feedback from the dust to the gas, since the gas 
density remains static. However, when the dust accumulates and 
grows forming a pile-up, the dust-to-gas ratio can reach values 
close to one. When the dust-to-gas density ratio is close to unity, 
the feedback from the dust to the gas is non negligible and the 
growth timescales can be faster than the drift timescales (Youdin 
& Goodman 2005), potentially leading to planetesimal forma- 
tion in the absence of self-gravity. Nevertheless, these streaming 
instabilities are powered by the relative drift velocities between 
dust particles and gas, therefore in the pressure maximum where 
the radial drift of the dust is reduced, these instabilities are un- 
likely and our assumptions are still valid. 

Figure 6 shows the dust-to-gas ratio for each considered case 
of the planet mass after 1 Myr of dust evolution. The maximum 
value that the dust-to-gas ratio is ~ {0.76, 0.36, 0.18) in the cases 
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of 1, 9 and 15 Mjup, respectively. This variation of dust to gas ra- 
tio with the planet mass can be explained considering that as the 
pressure bump gets wider, the dust accumulates and grows in an 
extended radial region. The dust-to-gas ratio is therefore higher 
for lower planet mass, when the pressure gradient is higher and 
positive in a narrower region (see Fig 3). In addition, we can no- 
tice in Fig. 6 that for the 15 Mjup mass planet, the dust-to-gas 
ratio reaches very low values at the location of the planet, for 
which the gas drag and drift is inefficient. Therefore the remain- 
ing dust in the inner region is just the initial dust that was initially 
in the inner part. If the initial size of the particles is taken smaller 
(e.g. 0.01 fiin), the sub-micron size particles would go through 
the gap and replenish the inner part. However, for the 1 and 9 
Mjup cases, the micron size dust from the outer region which are 
not trapped, can drift and be dragged by the gas and finally move 
through the gap. 




5RA { 



4. Observational predictions 

4.1. Emission maps 

In this section, we present continuum maps computed at a wave- 
length of 1.3 mm. We use the gas profiles described in Sect. 2.2 
and dust size distributions determined in Sect. 3.3 and shown in 
Fig. 4-middle panel and Fig. 5. We consider the same aspect ratio 
as in the FARGO and dust evolution simulations, and compute 
the volume density following: 



P(r, z) 



X exp 



IhirY 



(12) 



with Z, the dust surface density and h{r), the scale height. 
After retrieving the optical constants ' , we compute the opacities 
for silicate grains, following the Mie theory. The dust tempera- 
ture is then computed self-consistently after solving the radiative 
transfer with the RADMC3D code ^. 

We consider the cases of a = lO"-', a planet with masses of 
1 and 9 Mjup planets and orbiting at 20 AU, after 1 Myr. We use 
40° and 90°, as the disk inclination and position angle, respec- 
tively, and a distance of 140 pc. The images are convolved with 
a 0.10"xO. 10" beam as required for ALMA for typical sizes and 
distances of protoplanetary disks. As shown in Fig. 7, our model 
images are very similar to the transition disks observations in 
the (sub-)millimeter regime, that reveal a continuum emission 
as a resolved ring-like feature (e.g., Williams & Cieza 2011; 
Andrews et al. 201 la; Isella et al. 2012). In addition, after 1 Myr 
of dust evolution, we estimate a spectral index of ~2.7 by in- 
tegrating the flux over the entire disk at wavelengths of 1 mm 
and 3.0 mm. This value is in agreement with recent observations 
(Ricci et al. 2010). 

Due to the wide distribution of particles around the pressure 
bump (see, e.g.. Fig. 5), the brightness distribution is expected to 
change with wavelength. This feature is also reproduced by our 
model, as shown in Fig. 8. 

4.2. Tlie case of LkCa 15 

In this section, we apply and discuss our model in the context of 
the transition disk LkCal5, for which a wealth of observational 
data is published. Andrews et al. (2011a) presented a 870 fim 

' http://www.astro.uni-jena.de/Laboratory/Database/databases.html 
^ http://www.ita.uni-heidelberg.de/~dullemond/software/radmc-3d/index.shtml 



Fig. 8. Radial surface brightness in the 9 Mjup case, computed at 
160//m (crosses), 360/im (stars) and 870/im (diamonds). 



image obtained at high angular resolution (0.25") combining ob- 
servations from the Submillimeter Array (SMA) and Plateau de 
Bure Interferometer (PdBI). These observations of the dust con- 
tinuum emission confirms a ring-like structure at ~50 AU from 
the star They compared these observations with four possible 
models, concluding that their best fit corresponds to a gap pos- 
sibly opened by a 9 /Wjup companion at ~ 16 AU in a wedge 
disk, with a constant disk aspect ratio of /i/r = 0.05 and con- 
stant viscosity v = 10"^ (which corresponds to ff ~ 10"^). In 
Sect. 3, we demonstrated that, in the case of a 9 Mj^p planet and 
a - 10^^, the pressure maximum is located around ~ 1.4 rp 
from the planet (with r^, the planet location), which implies that 
the pile-up of dust (up to cm-sizes) would be located at around 
~ 38 AU. Therefore, to have an agreement with the depletion 
of dust opacities until ~ 50 AU, a bit more massive or/and fur- 
ther located companion is necessary to reproduce the ring-like 
structure as we will discuss below. 

Although not confirmed, Kraus & Ireland (201 1) claimed to 
have detected a ~ 6-15 Mj^p planet at 15.7 + 2.1 AU in LkCal5. 
Considering the upper limit for the mass and the location of the 
planet at 20 AU, the zero point for the pressure gradient would 
be located ~ 54 AU (Fig. 3). Hence, dust would accumulate at 
~ 54 AU distance from the star, in better agreement with the 
observed ring radius. 

Figure 9-top panel shows the vertically integrated dust den- 
sity distribution after 1 Myr of evolution in the case of a gap 
opened by a 15 /Wjup planet, embedded in a flared disk with 
a - 10^^ and the same parameters given in Tables 1 and 2. 
Considering this dust distribution. Fig. 9-bottom panel gives the 
1.3 mm continuum map, considering a disk inclination of 51°, 
a disk position angle of 64.4°, a distance of 140 pc, and con- 
volved with a beam of 0.21"x 0.19" (Isella et al. 2012). This 
image shares striking similarities with the ones published, with 
similar inner radii for the ring and similar surface brightness. 
We note however, that we do not reproduce the observed width 
of the ring, as after 1 Myr already, most of the outer disk par- 
ticles are retained in a very narrow region. Nevertheless, Pinilla 
et al. (2012) demonstrated that the presence of long-lived pres- 
sure bumps such as results from MRI can explain the presence 
of mm-size grains in the regions of the disk outside 50 AU. The 
combination of those pressure inhomogeneities and the interac- 
tion between a massive planet with the disk may explain the ra- 
dial extent. 
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Fig. 9. Vertically integrated dust density distribution at 1 Myr of 
evolution for a surface density with a gap created by 15 Mjup 
located in a fixed orbit at 20 AU in a flared disk and a - 10"^ 
(top-panel). Bottom-panel is the 1.3 mm continuum model map 
for LkCalS considering the dust density distribution shown in 
top-panel. The map is convolved with a beam of 0.21" x 0.19", 
corresponding to the maximum angular resolution achieved by 
Isella et al. (2012) Units are in Jy/beam. 

Unlike our models with an alpha-viscosity and a massive 
planet embedded in the disk, Isella et al. (2012) model a gas 
density profile which increases with radius for the inner part due 
to the extremely rapid increase of the viscous accretion velocity 
for smaller radius. They fit their CARMA 1.3 mm observations 
assuming no planet with the following gas surface density: 



r,\-y 



2(r) = 2;,(^) 



X exp < - 



1 



2(2 - r) 



2-r 

-I -1 

r, 



(13) 



where the viscosity is given by v(r) oc r''. Their best fit model 
corresponds to 7 = -2.15, r, = 60AUandE, ~ 10.37 g.cm"^, 
for a disk extending from 0.2-160 AU. Their gas density profile 
suggests that the pressure gradient would be positive from the 
inner radius until the characteristic radius r, and smoothly be- 
comes negative until the outer radius of the disk. This implies 
that the dust particles would have a positive drift velocity in the 



inner regions, and negative in the outer regions of the disk. They 
would therefore move outwards from the inner radius to r,, and 
from the outer radius inwards to r,, creating a wide concentration 
of dust around 60 AU. 

This complementary study, that uses a prescription of the 
surface density based on the viscous evolution of the disk, leads 
to the same conclusion as our model, i.e., a pressure maximum 
around 60 AU. In fact, the profiles are very similar in the range 
of radii in which the dusty ring is observed (see Fig. 1 in Isella 
etal. (2012)). 



5. Discussion 

Hydrodynamical simulations of planet/disk interactions com- 
bined with dust coagulation/fragmentation models can repro- 
duce the ring-like emission of some transition disks. 



5.1. Dust evolution 

To compute the evolution of dust, one has to account for four 
important physical mechanisms: coagulation, fragmentation, ra- 
dial drift and turbulent mixing. Coagulation is important for 
grains to grow to large sizes, as inferred from millimeter ob- 
servations. Fragmentation continuously replenishes the disk in 
small grains. If we do not consider fragmentation, in the region 
where the pressure gradient is positive, particles would accumu- 
late, grow above the fragmentation barrier and reach the maxi- 
mum size of ~2 m considered in our models. This would happen 
on timescales as short as ~ 0.1 Myr, and would result in a disk 
only made of large particles. As meter-size objects have St » 1, 
they would therefore not be coupled with the gas anymore, move 
at Keplerian velocity, and would not be trapped. These grains 
would in addition have very low opacities, which would result in 
the underestimation of millimeter fluxes in the millimeter 

Radial drift is essential to distribute the grains over the disk 
extent. Birnstiel et al. (2010b) showed that if the radial drift is 
neglected, dust grains could acquire millimeter sizes, but would 
lead to over-predictions of mm-fluxes. In our model of LkCal5, 
we overcome this issue and obtain flux levels similar to the ones 
measured by Isella et al. (2012). 

Finally, turbulent mixing provides the necessary relative ve- 
locities for grains to fragment, directly influencing the maximum 
particle size (Eq. 4), and distribute the grains of various sizes all 
over the disk. In Sect. 3.3, we showed that the efficiency of par- 
ticle trapping depends on the viscosity of the disk. For a shallow 
gap (see Fig 4), we concluded that when the pressure gradient is 
not high enough, drag forces and diffusion do not allow particles 
to be trapped. 

Our model can however still be refined. Recent collision ex- 
periments with silicates show that particles should also bounce 
as an intermediate step between sticking and fragmentation 
(Guttler et al. 2010). In fact, Windmark et al. (2012) showed that 
with the insertion of cm-sized seeds, considering mass transfer 
processes and bouncing, at a distance of few AU, small dust can 
grow to ~ 100 meters on timescales of 1 Myr A similar process 
could happen to mm-sized particles in the outer regions of the 
disks. However, in the presence of ices, the need to account for 
bouncing effects is still highly debated (e.g. Wada et al. 2011), 
and we do not consider it in this work. 
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5.2. Different gap extents in tine dust and tine gas 

The presence of a planetary or sub-stellar companion has been 
suggested in many studies to explain the properties of transition 
disks. However, studies based on hydrodynamical simulation of 
gap opening in a gaseous disk have shown that a single planet 
could not open the large gap inferred from spatially resolved 
millimeter observations (e.g. Dodson-Robinson & Salyk 2011; 
Zhu et al. 2011). Most studies in fact consider a maximum ra- 
dius for the gap of 5 Hill radii. According to the definition of 
the Hill radius, the increase of the planet mass does not pro- 
duce a significant widening of the gap, and none of these models 
with a single planet can reproduce the observed characteristics 
of transition disks. In this study, we include for the first time 
the modelization of dust growth and evolution, in addition to the 
hydrodynamical simulations. We predict grain size distributions 
(and therefore brightness distributions) and show that the gap 
edges are different in dust and gas, and that therefore the extent 
of the gap due to a single planet is larger than expected based 
on gas simulations. We demonstrate in Sect. 3 that for a given 
planet, dust particles are trapped further out than the location of 
the edge of the gap in the gas. For example, in the specific case 
of a 15 Mjiip planet at 20 AU, the outer edge for the gas would 
be located at 37 AU, while the dust ring would be at 54 AU. 

For Mp > 5 Mjup, the density waves produced by the planet 
are still changing with time after few thousands orbits. Kley & 
Dirksen (2006) showed that the disk becomes eccentric with the 
presence of a planet with a mass Mp > 3Mjup (for a constant 
viscosity given by v = 10"^), which may cause a change in the 
equilibrium state. However, the disk precesses slowly enough 
that the eccentric pattern appears to be nearly stationary, and for 
a critical mass of 5 Mj^p, the eccentricity of the disk reaches a 
threshold. In fact, we demonstrated in Sect. 3.2 that for 1 and 
3 Mjup the dust would accumulate at 7 rn from the location of 
the planet r^, and at 10 rn from rp for 5, 9 and 15 Mjup. We 
attribute this difference in the locations of the zero point of the 
pressure gradient to the fact that the disk becomes eccentric for 
planets whose mass is higher than 3 Mjup. As a result, we have 
the pressure maximum further away for more massive planets, 
when the disk is eccentric, than for less massive planets. 

5.3. The inner disk 

Near-infrared excesses are measured in some transition disks 
(sometimes called pre-transition disks), and indicate the pres- 
ence of dust close to the star, that can be spatially resolved in 
a few cases (e.g., Olofsson et al. 2011; Benisty et al. 2010; Pott 
et al. 2010). Besides, these objects show signs of on-going accre- 
tion (see e.g Espaillat et al. 2007). In this paper, we do not treat 
the inner disk because our prescription for it is not accurate. In 
order to have self-consistency between gas hydrodynamical sim- 
ulations and dust evolution models, we considered for the kine- 
matic viscosity v = o- c, h. With this viscosity and the two con- 
sidered assumptions for the gas density profiles, the accretion 
rate M = 3 tt v Z is independent of the radius. With the disk 
and stellar parameters considered for the two cases (Tables 1 and 
2), M is of the order of ~ 10"'' Mq yr"', about an order of mag- 
nitude lower than the measured values in some transitional disks 
(~ 10-** Mo yr-i; Calvet et al. (2005); Espaillat et al. (2007)). 
The low model value implies that even considering open bound- 
ary conditions, the gas surface density, in both inner and outer 
regions, only slightly decreases over the 1000 orbits (Fig. 1). In 
addition, the insertion of a very massive planet in the disk leads 
to strong density waves that may lead to unrealistic effects at the 



inner boundary (e.g., Crida & Morbidelli 2007). Therefore, our 
prescription for the inner disk is not accurate. 

Besides, it is important to notice that inside the ice-line 
the absence of ices can change our results. In fact, accord- 
ing to laboratory experiments, the fragmentation velocity varies 
with the material properties (e.g. Gundlach et al. 2011). In the 
presence of ices, the fragmentation velocity reaches a value of 
Vf ~ 1000cm s"', the value that we consider for the outer disk. 
However, this threshold decreases in the inner AUs, inside the 
ice-line. This directly impacts the maximum particle size (Eq. 
4). 

Finally, the radial drift in our model leads to empty the inner 
disk in timescales that can be shorter than the age of some tran- 
sition disk with NIR excess (e.g., HD100546). Replenishment 
from the outer disk with small grains passing through the gap, or 
from planetesimal collisions (Krijt & Dominik 201 1) is needed. 

6. Summary and conclusions 

The sub-Keplerian radial velocity of the gas in protoplanetary 
disks makes millimeter-size particles in the outer regions of the 
disk exposed to a rapid inward drift, implying that they mi- 
grate towards the star on timescales shorter than a million years. 
Therefore, in planet formation, rapid inward drift is one of the 
main issues with models to form planetesimals. The idea of the 
presence of pressure bumps in protoplanetary disks has been pro- 
posed as a solution to stop the rapid inward drift. With the pres- 
ence of a massive planet in a disk, a pressure bump is created 
in the outer edge of the cleared gap. Particles may experience a 
positive pressure gradient and as a result, do not drift anymore 
and not experience the high-velocity collisions due to relative 
radial and azimuthal drift. Nevertheless, particles can still frag- 
ment due to turbulence motion, and the resulting micron-size 
particles are less easily trapped and they can finally drift inwards. 

Some transition disks reveal gaps that can result from the 
presence of a massive planet, making them potentially interest- 
ing laboratories for studying dust growth under the favorable 
circumstances of having a significant pressure bump. In this pa- 
per, we combine two-dimensional hydrodynamical simulations 
for the gas with one-dimensional coagulation/fragmentation dust 
evolution models, to study how dust evolves in a disk which gas 
density has been disturbed by a massive planet. 

We investigate the influence of the disk geometry, turbulence 
and planet mass on the potential trapping of particles. The disk 
geometry does not have a significant effect on the gap opening 
process. For a 1 Mjup planet, there is an important influence of 
the a turbulence parameter on the depth of the gap and has im- 
portant consequences for the trapping of particles. Unlike the 
case of Q- = 10"^, with a - 10^^, the particles are trapped at a 
distance of ~ 0.5 from the planet position rp. While the gas 
gap has a radius of ~ 0.35 rp (or 5 Hill radii), the dust is retained 
at a further distance. Without considering turbulence mixing for 
the dust dynamics would produce trapping of particles for cer- 
tain sizes regardless of a value provided that the pressure gra- 
dient is positive. In the case of 9 Mjup planet, the gap depth is 
independent of the turbulence and as a consequence the pressure 
gradient behaves similarly for all a values. The dust particles are 
trapped at a distance of 1 .4 rp from the planet orbit rp. The planet 
mass strongly influences the location where the dust is retained 
and the width of the dust bump. 

Observations at millimeter wavelengths reveal some transi- 
tion disks with very wide gaps. We show that the location where 
the dust piles up does not coincide with the gap outer edge in the 
gaseous disk. This mismatch suggests that multiple planets may 
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not necessarily be needed to account for the observed large opac- 
ity holes. We reproduce the main observed features of the disk 
around LkCalS, in which a companion was recently detected. 
Reproducing asymmetric features in the disk, such as the ones 
found in HD135344B (Brown et al. 2009) is subject of future 
work, as they will be well detectable with ALMA. 
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